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Abstract 

We perform classical molecular dynamics simulations to demonstrate the application of kinked 
silicon nanowires (KSiNWs) as nanoscale springs. The spring-like oscillation in gigahertz frequency 
range is successfully actuated using a similar procedure as the actuation of a classical mass spring 
oscillator. We detect the spring-like mechanical oscillation and some other low-frequency oscilla- 
tions by the energy spectrum analysis, where a dimensional crossover phenomenon is observed for 
the transverse mode in KSiNWs with decreasing aspect ratio. Our findings shed light on the elastic 
properties of the KSiNW and open a way for its application in nanomechanical devices. 
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The kinked silicon nanowire (KSiNW) was synthesized by Tian et al. in Lieber's group 
at Harvard University in 2009 1 -. Owing to its peculiar kinking structure, this nano struc- 
ture is viewed as one of the promising building blocks in bottom-up integration of active 
devices. It is able to manipulate the arm length of the kink by controlling the growth time. 
The change of growth direction at kinks is the characteristic of KSiNWs; thus great exper- 
imental efforts have been devoted to investigating the growth mechanism of the kinks in 
silicon nanowires 2 - - -. Besides these experimental works, Schwarz and Tersoff have proposed 
a theoretical model to interpret the growth mechanism of the KSiNW- 1 ^. In their model, 
the kink formation is attributed to the interplay of three basic processes: facet growth, 
droplet statics, and the introduction of additional facets. Stimulated by KSiNWs, several 
groups have synthesized kinks in other nanowires, such as In 2 03 multikinked nanowires^, 
kinked germanium nanowires 1 ^, germaniumsilicon axial heterostructure with kinks 1 ^, and 
kinked ZnO nanowires 1 ^. As long as significant achievement has been gained for the kink 
growth mechanics, a meaningful task would be the investigation of possible practical appli- 
cations of such material. A recent experiment has demonstrated the application of KSiNWs 
as nanoelectronic bioprobes.— Hereby, we demonstrate its utility as a nanoscale spring in 
mechanical device. 

In this paper, we study the mechanical oscillation of KSiNWs of various aspect ratio by 
molecular dynamics (MD) simulation. We demonstrate a successful actuation of the spring- 
like oscillation in the KSiNW by mimicking the actuation procedure for a classical mass 
spring system. We perform an energy spectrum analysis for the mechanical oscillations, 
where the spring-like oscillation and some other low-frequency oscillation modes are suc- 
cessfully assigned to corresponding resonant peaks in the spectrum. From energy spectrum 
analysis, we also find the dimensional crossover behavior for the transverse oscillation mode 
in KSiNWs with decreasing aspect ratio. 

Fig. H](b), (c), and (d) show the configuration of three samples in present study: sample 
1 (SI), sample 2 (S2), and sample 3 (S3). In all samples, the nanowire growth direction 
changes from one < 100 > direction to another < 100 > lattice direction, resulting in a 
kinking angle of 90°. All of the four side surfaces are the {100} lattice surface, enclosing a 
square cross section. The cross section is of 27 x 27 nm 2 . For SI, S2, and S3, the aspect 
ratio (length to thickness ratio) is 3.6, 3.0, and 2.4; while the number of silicon atoms is 
8 x 10 4 , 6 x 10 4 , and 4 x 10 4 . 
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All MD simulations were performed using the publicly available simulation code 
LAMMPS^, while the OVITO package was used for visualization^. The interaction between 
silicon atoms is described by the Stillinger- Weber potential^. Owning to its efficiency and 
accuracy, this empirical potential has gained a wide application in the simulation of silicon 
and other similar valence bonded systems. The Newton equations of motion are integrated 
using the velocity Verlet algorithm with a time step of 1 fs. In the growth direction, the 
left end (black online) is fixed in all simulations, while the right end is free. Free boundary 
condition is applied in lateral directions. 

To actuate the mechanical oscillation of KSiNWs, we first learn a lesson from the classical 
approach to generate the oscillation of a harmonic mass spring oscillator shown in Fig. [T] (a). 
The point mass on the right end is displaced away from its equilibrium position Z by an 
external force F to an updated position I. The system starts to oscillate harmonically after 
releasing F. The oscillation amplitude is A — I — Iq. The maximum strain inside the spring 
during the oscillation is e m = A/ lo- 
in a similar way, we actuate the mechanical oscillation of the KSiNW in following four 
steps. Firstly, the structure is relaxed to the minimum energy state using the conjugate gra- 
dient algorithm. Secondly, the system is stretched gradually by adding a ramp displacement 
that scaled linearly from zero at the fixed left end to a maximum value at the free right 
end. The displacement incremental results in a strain rate of e = 10 9 s _1 , which is a typical 
value in MD simulations^. This displacement incremental procedure is stopped after the 
tension inside the system reaches the required value, e m . In following simulations, e m = 0.01 
or 0.02, which are small enough to avoid plastic deformation in the KSiNW. Thirdly, the 
right tip (blue online) in the tensile system is clamped, and the system is thermalized by 
the Nose-Hoover— thermostat to a constant temperature of 4.2 K within the NVT (i.e. 
the particles number N, the volume V, and the temperature T of the system are constant) 
ensemble for 20 ps. We focus on the liquid helium temperature because this temperature is 
commonly utilized in experiments involving nanoscale mechanical oscillations (eg. 90 mK 



or 4.0 K in Ref. 



2ll ). Fourthly, the right tip is set free and the system is allowed to oscillate 



freely within the NVE (i.e. the particles number N, the volume V, and the energy E of the 
system are constant) ensemble for 4 ns. The trajectory and energy time history from the 
final step are then used in the analysis of the mechanical oscillation for the KSiNW. 

Fig. [2] shows the energy time history in KSiNW SI during the simulation within NVE 
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ensemble, with e m = 0.01. All energies are with reference to their value at t = 0, i.e at 
the very beginning of the step 4 in the above actuation procedure. The kinetic (potential) 
energy is above (bellow) the x axis; while the total energy locates exactly at the x axis. Cor- 
responding to the mechanical oscillation, the energy oscillates between kinetic and potential 
energy, and an excellent energy conservation has been achieved. This nice energy oscillation 
behavior exhibits a good spring-like behavior of the KSiNW. We note that the frequency of 
the energy oscillation (/ = 31.26 GHz) is twice of the mechanical oscillation frequency. Two 
cosine functions (red online) are also shown as guides to the eye. 

For quantitative analysis of the spring-like oscillation, we record the vibrational displace- 
ment of the right tip u(t) = (u x ,u y ,u z ) = r(t) tip — r^ ip . u is averaged over all silicon atoms 
within the tip on the right end. These vibrational displacements are plotted in Fig. [3] (a), 
(b), and (c) for the three studied samples SI, S2, and S3. The actuation strain here is 
e m = 0.01. In all samples, a steady oscillation in x direction has been successfully actuated. 
This oscillation reflects the spring-like behavior of the KSiNW. The oscillation in z direction 
only induces some neglectable fluctuations. However, we can observe some obvious oscilla- 
tions in the y direction, although our original goal is to actuate oscillation solely in u x by 
displacing atoms in the x direction. This is due to the so-called 'artificial' effect as discussed 



in Ref. 
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In that work, it was found that lots of other oscillation modes will be excited 
simultaneously if the mechanical oscillation is not actuated properly according to its funda- 
mental vibration modes. In present work, the mechanical oscillation is actuated by a linearly 
ramp displacement. Due to surface and finite size effects, the linear distribution shape does 
not correspond to the fundamental vibration mode in the KSiNW, so the artificial effect 
leads to the excitation of lots of other oscillations besides the spring-like mechanical oscilla- 
tion. From the u x curve in Fig. [3] (a), the mechanical oscillation frequency is extracted to be 
/ = 15.63 GHz, which is exactly half of the frequency in the energy oscillation in Fig. [2j It 
indicates that the spring-like oscillation in x direction makes most important contribution 
to the oscillation in the kinetic/potential energy, even though its oscillation amplitude is 
somehow smaller than the oscillation in u y . We point out some roughness on both u x and 
u y curves, which are actually due to higher-frequency oscillation modes. 

From Fig. [3], we have analyzed the spring-like oscillation based on the vibrational dis- 
placement. In this analysis, we have encountered some disruption induced by lots of other 
oscillation modes resulting from the artificial effect. The artificial effect can be eliminated 
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by actuating the mechanical oscillation complying with the fundamental vibration mode. 
However, it is impossible to get the fundamental vibration mode for those KSiNWs here, 
because their lattice dynamics matrix can't be solved due to huge degrees of freedom. On 
the other hand, the oscillation in kinetic/potential energy shown in Fig. [2] disclose nicely the 
spring-like mechanical oscillation. Further more, a direct analysis from the vibrational dis- 
placement is much more difficult, even though some significant progress has been achieved^; 
while it has been a standard setup in the experiment to analyze the mechanical oscillation 
from the energy spectrum. To this end, we should also study the spring-like mechanical 
oscillation in KSiNWs from their energy spectrum. 

Fig. H] shows the energy spectrum for all samples, which are actually the Fourier transform 
of the kinetic energy time history—. Each peak in the figure corresponds to a special oscil- 
lation mode in the system. We note that the oscillation mode here is essentially equivalent 
to the concept of phonon vibration mode in the lattice dynamics theory. There are three 
low-frequency branches for vibration modes in KSiNWs: one for longitudinal acoustic (LA) 
mode and two for transverse acoustic (TA) modes. In KSiNW SI, a strong peak locating at 
/ = 31.26 GHz is assigned to the first lowest-frequency LA (LAI) mode. On the left side of 
this LAI peak is a smaller peak at f — 18.60 GHz corresponding to the first lowest-frequency 
transverse mode in y direction (TA y l). Both LAI and TA y l peaks are successfully assigned 
by comparing the position of these peaks to the frequency extracted from Fig. |3j We have 
also denoted the two higher-frequency modes in x and y directions, i.e LA3 and TA J/ 2, at 
/ = 88.69 and 84.16 GHz. The frequency of LA3 is threefold of the frequency for LAI, 
because the frequency of the LA mode is linearly proportional to the wavelength. Hence the 
frequency is linearly proportional to its mode index. There are only LA modes with odd 
index in such quasi-one-dimensional structure with one end fixed and the other end free^. 
The frequency of TA y 2 mode is about fourfold of that for the TA^l mode. It is because 
of the flexural nature of the transverse modes in one-dimensional rod-like structure of high 
aspect ratic^ir— , where the frequency is proportional to the square of the mode index. 

For KSiNW S2, we also successfully assign the three low-frequency modes, LAI, LA3 
and TAyl, to the three strongest peaks at / = 53.17, 156.40, and 32.90 GHz, respectively. 
However, we do not observe a peak corresponding to the TA y 2 mode, which should occur 
around / = 131.60 GHz. Similar as S2, we can assign the three low-frequency modes, LAI, 
LA3 and TA^l, to the three strongest peaks in KSiNW S3 at / = 117.78, 359.54, and 
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82.49 GHz. Similarly, we do not see a peak corresponding to TA y 2 mode, which should 
locate around / = 329.96 GHz due to the flexural nature of the transverse mode. However, 
there is an obvious peak at / = 182.63 GHz (denoted by a star). Inset of the figure shows the 
Fourier transform for the vibrational displacements u x , u y , and u z . There is a strong peak 
at / = 92.51 GHz for u y , indicating that the resonant peak at / = 182.63 GHz in the energy 
spectrum should be assigned to a particular oscillation mode in the y direction. It means that 
this peak corresponds to the TA y 2 mode in S3, which is contrary to the predicted position 
/ = 329.96 GHz from its one- dimensional flexural nature. This contrary actually reflects 
the three-dimensional bulk property of the KSiNW S3. As well known, in bulk system, the 
transverse mode is not flexural anymore. Instead, its frequency is also linear proportional to 
the wavelength, so its frequency is also proportional to the mode index number. According 
to this bulk prediction, the frequency of the TA y 2 mode should be twice of the TA y l mode, 
i.e / ~ 164.98 GHz. This value is very close to the peak at / = 182.63 GHz in the energy 
spectrum. This finding discovers the dimensional crossover (from one-dimensional to three- 
dimensional) for the transverse mode in KSiNWs with decreasing aspect ratio. 

To detect possible effect from the actuation strain e m , we have also performed simulations 
using e m = 0.02. The obtained frequencies of the spring-like mechanical oscillation are 
compared with that from e m = 0.01 in Tab. HI There is only slightly enhancement of the 
frequency by using larger actuation strain. This enhancement effect is weak because of the 
free boundary condition on the right end, which leads to the absence of the effective strain. 
It is better to fix both ends to induce effective strain for applications where strong frequency 
enhancement is desirable.— 

To summarize, massive MD simulations have been performed to investigate the spring- 
like mechanical oscillation of KSiNWs with three different aspect ratio. The spring-like 
oscillation is actuated by the method borrowed from the actuation of a classical mass spring 
oscillator. The energy spectrum is utilized to analyze the mechanical oscillation of KSiNWs, 
where the spring-like behavior is detected, accompanying with some other low-frequency 
oscillation modes. A dimensional crossover is observed for TA y mode in KSiNW samples 
with aspect ratio decreasing from 3.6 to 2.4. 
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FIG. 1: Configuration of the three KSiNW samples in present work, (a) The equilibrium position 
of a spring connected to a point mass. The right end is stretched by an external force F (with 
fixing left end). This mass spring oscillator starts to oscillate after F is removed. From (b) to (d): 
the three studied KSiNW samples SI, S2 and S3. Left end is fixed. Tips on the right ends are in 
blue online. Cartesian coordinate system: x is along axial direction from left to right, and z is the 
viewing direction of structures shown here. 
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MD time (ns) 

FIG. 2: Energy oscillation due to the mechanical oscillation in KSiNW SI. Energies are with 
reference to the value just before the actuation of the mechanical oscillation (i.e t=0). The curve 
above (bellow) zero is the kinetic (potential) energy. The energy oscillates between kinetic and 
potential energy at frequency / = 31.26 GHz, corresponding to the mechanical oscillation. The 
total energy is the horizontal line at zero, demonstrating a good energy conservation in the NVE 
ensemble. Two cosine functions (red online) are used as guides to the eye. We note that the energy 
oscillation frequency is double of the mechanical oscillation frequency. 
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FIG. 3: Vibrational displacements of the tip at the right end in KSiNWs. Top to bottom: results 
for SI, S2, and S3. u x , u y , and u z are the vibration in three Cartesian directions. 
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FIG. 4: The energy spectrum, i.e the Fourier transform of the kinetic energy time history. The 
first several low- frequency oscillation modes are denoted for all of the three studied samples. With 
decreasing aspect ratio from SI, S2 to S3, the second transverse mode in y (TA y 2) undergoes a 
transition from one-dimensional nanowire behavior into three-dimensional bulk behavior. Inset: 
Fourier transform of u x , u y , and u z for S3. 
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